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ABSTHACT 


The  scattering  of  a plane  ware  by  a spherical  bubble  immersed 
in  a homogeneous  perfect  fluid  of  infinite  extent  has  been  examined. 
Expressions  for  (l)  the  average  sound  intensity  of  the  scattered  wave,  and 
(2)  the  average  total  intensity  of  the  entire  sound  field  have  been  obtained. 
Numerical  calculations  have  been  carried  out  for  the  case  of  long  wave- 
lengths and  for  distances  close  to  the  bubble.  Several  curves  have  been 
plotted. 
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A STUDY  OF  THE  SCATTERING  OF  SPUED  WAVES  IN  FLUIDS 
PARTl  - SCATTTOIHQ  OF  A PUBS  EA7g  BT  k SPHERICAL  BOBBLE 

I.  IHTRODUCTIOM 

The  scattering  of  a plans  wave  by  a spherical  bubble  immersed  in 
a homogeneous  perfect  fluid  of  infinite  extent  has  been  examined.  An  ex- 
pression for  the  average  sound  intensity  of  the  scattered  wave  very  close 
to  the  bubble  has  been  obtained.  The  important  case  of  long  wavelengths 
or  snail  bubble  radii  is  of  particular  interest  and  has  been  given  special 
attention.  numerical  results  are  presented  in  the  fora  of  polar  diagrams 
for  several  different  values  of  r/ec  and  Zt>k/K  , where  r is  the  distance 
froa  the  center  of  the  bubble,  a.  is  the  bubble  radius  and  A is  the  wave- 
length. 

In  addition  to  the  average  intensity  of  the  scattered  wave,  an 
expression  for  the  average  total  intensity  of  the  entire  eourd  field  has 

( ') 

been  obtained.  This  result  is  cohered  with  the  average  sound  intensity 
of  a plane  wave  travelling  in  an  infinite  homogeneous  fluid  with  no  bubble 
preeent.  It  is  found  that  the  influence  of  the  bubble  on  the  average  total 
sound  intensity  is  negligible  when  Zrra/K.  is  snail  and  p/ft.  »I0.  Thus, 
the  average  total  sound  intensity  differs  very  little  froa  that  due  to  a 

j 

plane  wave  travelling  in  an  infinite  homogeneous  fluid  with  no  bubble 
present.  For  snail  values  of  r/a. , however,  the  press noe  of  the  bubble 
seriously  alters  the  magnitude  of  the  average  total  sound  intensity  froa 
that  due  to  a plane  wave.  Polar  diagrams  of  the  ratio  of  the  average 
total  sound  intensity  to  the  average  intensity  of  the  incident  plane  wave 
are  presented  for  several  different  values  of  r/a. . 
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II.  DISCUSSION 


Consider  a plane  wave  incident  upon  a spherical  babble  of  radius 
A.  iaaereed  lp  an  infinite  homogeneous  perfect  fluid  characterised  by 
density  f>  and  velocity  C . It  is  desired  to  find  expressions  for  the 
average  sound  intensity  of  the  scattered  wave  and  the  average  total 
intensity  of  the  entire  sound  field  very  close  to  the  bubble.  The  special 
case  when  the  ratio  of  wavelength  to  sphere  radius  is  large  is  to  be  given 
special  attention.  Large  values  of  this  ratio  laply  that  either  the  wave- 
length is  long  or  the  sphere  radius  is  wall. 

The  pressure  equation  for  a sound  wave  is  given  by 


ft!  ’ 


(1) 


where  p is  the  excess  pressure  and  C is  the  velocity  of  sound.  If  the 
pressure  dependence  on  tiae  is  of  the  fbrw  exp(-cwt),  then 


v>+t>=°  > 

where 

i-o  2nV  2W 

The  scalar  wave  equation  in  spherical  coordinates  is  given  by 


(2) 


(3) 
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wher*  r*  t $ and  are  shown  in  Figure  1,  The  problem  under  discussion 
possesses  symmetry  about  the  polar  axis  so  the  incident  and  scattered  wanes 
are  independent  of  ^ . In  order  to  solve  the  problem  it  is  necessary 
to  find  solutions  of  the  scalar  wave  aquation  that  represent  an  incoming 
plane  wave  and  an  outgoing  scattered  wave.  In  addition  to  those  require- 
ments, it  is  necessary  that  the  solutions  satisfy  the  proper  boundary 
conditions  on  ths  surfscs  of  ths  bubble. 

The  expression  for  a plane  wave  travelling  to  the  right  along  the 
polar  axis  in  teres  of  spherical  waves  is  given  by 

pfm  i.  «• 

or  (5) 

ft-  r.Z  (z 

*H»0 


where  are  Legendre  Fur.ct  .ons  of  the  First  lind  and  ^Ck*')  are 

Spherical  Bessel  Functions.  The  expression  for  the  wave  scattered  by  a 
sphere  of  radius  A.  whose  center  is  at  the  polar  origin  is 


(6) 


where  njkr)  are  Spherical  Neumann  Functions.  The  combination 
[fm(tr)+i A.M]  represents  an  outgoing  wave. 

The  appropriate  boundary  condition  for  the  bubble  is  that  the 
total  pressure  on  its  surface  be  tero.  Therefore,  the  condition  to  be 
satisfied  when  is  that  the  sum  of  ths  incident  and  scattered  pres- 


sures be  ssro 


-I 'Mi' 


P,+f>r=°  ■ O' 

The  constants  C.  can  be  found  by  substituting  equations  (5)  and 


(6)  Into  equation  (?)<  Substitution  yields 


(?*"+*)  i"  (<■**&)  £m(k*)  t 


- uot 


ZC  fa)] e' 


- c out 


mO 


f?  6-  +l)i  ^ (U)  +C„  [fl  (U)+cn„  U*)]  • 0. 
[ ft*)  + Ur  )J  ^ - l (kr)  t * ^ ) ) 


where 


AJkr)  -j [}„»•)? +[»-0*)Y 

^(tr)  = 4U' Mtl  . 

- »~f  (kr) 

By  asking  use  of  the  relations 

fM  (kr)  = /f*  (kr) 

r»H  _ C?6»+t) 

L — C a 
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the  constant*  C,~  can  be  written 


n p/  \ -1' . v 

C,H  ~ ~~  *o  (2>n  + V C oc„  (Ihk)  . 

The  expression  for  the  scattered  pressure  ware  becomes 


h~-PaZ (** ') e i[*-(u) 

r*  - G 


•A  (lr)  *+■*  0<m  (ix) 


- iuit 


(ill) 


(15) 


by 


The  radial  coc^>onent  of  Telocity  for  the  scattered  ware  is  given 


/ 


2El  . 


- wp  2^ 


(16) 


By  substituting  equation  (6)  into  (16)  the  expression  for  the  Telocity 


of  the  scattered  waee  be coses 


21  ^ cr,~ 1 Lt  ^ 


• Icjt 


or 


oo 


^ = <>rZ  ^ f? Y«m»)[frfr)  ±i»'e  e^]pcut 

It  is  convenient  to  introduce  the  following  notation! 


(17) 


(18) 


* 


-6- 


By  using  the  appropriate  recurrence  for  aulas  it  can  be  sheen  that 

4 =rho I frv*  ”>4 1 ^ ^ ~P2t~ ' 


(19) 


and 


a (kr)  - 1*£  lir)~t  in  M 

l>*p-,(kr)-(p+i)n^'(kr) 


(20) 


Substituting  (lli)  and  (18)  into  (17),  the  expression  for  the  radial  Telocity 
of  the  scattered  wave  beooaes 

*«.  - - 4-  S ft,„j  e-  ‘hM  <M  -*  *> 

/#C 

• i?,  ft-)  ^ «>  (U)  f?(o*9 Je~  ^ . 


(21) 


In  a similar  way  it  is  found  that  the  radial  component  of  Telocity 
for  the  incident  plane  ware  is 


where 


7t : ^ 


-iwt 

e e 


(22) 


(23) 
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The  angular  component  of  Telocity  for  the  scattered  vara  la  given 


by 


/ 


ULf  - — • 

v iwf  r 

Differentiating  equation  (15)  with  reaped  to  & and  oombining  teims 

c*») *" ' [^(U)~  ^(tr)  * ***)J 


(2M 


o 


ft1*)  «>  ft*)  ^ ^ c*e  6?)  c 


-tut- 


(25) 


where 


The  angular  oo^onent  of  Telocity  for  the  incident  plane  wave  la 

o« 


fa-o 


i%Q>+')_Cufc 

t 


(26) 


The  average  in  tensity  of  the  scattered  wave  can  be  found  by 
■ultiplying  the  real  parts  of  the  expressions  for  the  scattered  pressure 
and  scattered  radial  Telocity  and  aTer aging  over  ties.  Note  that  the 
average  scattered  intensity  is  the  average  rate  at  which  energy  is  trans- 
mitted per  square  centimeter  in  a radial  outward  direction.  The  real 
parts  of  the  expressions  for  ft  and  Us  are  given  by 


. > 
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Therefore,  the  intensity  of  the  spatt^ed  ware  bee  cnees 
ffjL  (“*)  » -j~-  + 

6 fiat 

( cn&)f£  (cntyA^  (tb)  Bf>(tr>-)  O+[*^0*)  - - £*i  -f  ui  J 

p)  (Ib)--ffa/)  -h  * 

(29) 

By  using  the  appropriate  trigonometric  formulae  and  omitting  all  terms 
whose  time  average  is  tero, 

Co*-/**,  (U)~  £ » + “*]  *++[«/£*)  ~fy  Or)-  fO'+/)±u><  ] 

- f <^[**(^-*/k)-«Jlr*)+fy(kb)-i>*+'%(f>rl)]  (30) 

Substituting  the  above  result  into  equation  (29)  it  is  found  that  the 
average  intensity  of  the  scattered  wave  is 


t ■ — 

2?  ~ ( U)  tu+Letp  (h) 

>.->  r o p - 0 

• fcfaa)  tffcj*  o)A»  (kr)Bt  (lr ) 

■CJ+[*~(U)  (U) -«m  (kr)  -hfy  d>) -J  + T ty* ')]  , 


(31) 


where  1C0  w Zfa  Is  the  average  intensity  of  the  incident  plane  ware. 
Vote  that  if  c.g.s.  units  are  used,  the  intensity  is  Measured  in  ergs 
per  second  per  square  centimeter. 

The  average  total  intensity  of  the  entire  sound  field  can  be  found 
by  taking  the  vector  sun  of  the  average  total  radial  intensity  and  the 
average  total  angular  intensity.  The  average  total  radial  intensity  can 
be  found  by  Multiplying  the  real  parts  of  the  expressions  for  the  total 
pressure  and  total  radial  velocity  and  averaging  over  tine.  Note  that 
a positive  value  of  the  average  total  radial  intensity  is  the  average 
rate  at  which  energy  is  trans Kitted  per  unit  area  in  a radial  outward 
direction.  The  average  total  angular  intensity  can  be  found  by  multiply- 
ing the  real  parts  of  the  expressions  for  the  total  pressure  and  total 
angular  velocity  and  averaging  over  time.  In  this  case  a positive  value 
of  the  total  angular  intensity  is  the  average  rate  at  which  energy  is 
transmitted  per  unit  area  in  the  direction  of  increasing  $ . 
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The  total  pressure,  to«al  radial  velocity  and  total  angular 
velocity  ere  given  by  the  following  expressions: 

t*tmh 

~ K y~ C2,r'+l')Rt(r*+-$) ^ (kr)  e ~ € 


nf  *0 


^ y (2*'+ 1)  .a^  ( kr)  (k<) 

— * . , 

~ <■'  [*„  (k*)-  (kr)  - r ” J <T  CU)^ 

* d 

*th  = 

P.  ■?/-  . A n,  ^ \ ;?(/•*')  -iwt 


-~Zr2Af')^eWWe‘ 

• bs  o " 

t 


>=o 


-Lfa(U)-Pfdr)-  -iwt 

• e,  r e • 

t/L  ty*)  %(“**■ 9)/y^c 


(32) 


(33) 


■f'H  .«!,<■ 

e 


/e=° 


(<**&) A r MA 

/°C  />so 

- L [k  . (k^—^Ckr)  - y((  ^ J -tu/t 


(k<x) 


(3U) 
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The  real  parts  of  these  quantities  are  glren  by 

•C 

m^O 

_ p *y~ (2*+i) cb) 


M * 0 

• urz,  [cX*  (k*)  — <*-.  ( kr)  - y *,  + 

p.m  ^ 

h-o 


(35) 


^(U*r)  " ~~pZ  / fy* tiff (*+*8) £t(b)  f+f-ffftfy+wt] 


_ o* 

— ^ (cs*. 9 )iy  (kt)4**^0<p  CL^) 

fro 

•c Q+* [<Xp(ko)—  fy(lr)  — 4-i) 

P / x° 

fro 


(36) 


4 £ X (2^'(a*»)At(b)  «t  a*.) 

r=o 


(37) 
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The  total  radial  intensity  beconea 

» (76.  (fa )/6L  (ut,) 

^ rr%0  f>- O 

• /_  +U)  t J UH-f  l)  f ookj 

i oq  »o 

' »rr:  0 p=  0 

. >tA.  <xM //aJ/^vv  <*,  A)  CA+[ cx„  (la)-** (lr)  - 2*^  tut] 

. cto-y-  fa(U)  -ft(tr)-$(^')+i*t] 

x »■  CM 

- ~ /~ ^ fa+i)(2f+*)FZ(<*(ji\l^(Q*b)jm  M&[(kr)di*.tXf(i+) 

*• O f*D 


• C++[-z>n  + ojt]  (^[of/U)  ~ft(lr)-j('l»  + l)±^] 
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('  ) 


^ toUl  angular  intensity  bacons 


f pL  ~ 

~ ^ h (2~-fl](2^Z^e)P/Ue)t  a.  )^(i, ) 


F0X  / vT 

’^/T4  p-c  (2~^+‘^ ^)%Yc*8)L  (lr)At  (i,) 

a*)  c^[<k,  (U)-*.,  a.)  - p +ut] 

■ CM  /V^)-~, rt-j  ,) v-^t ] 

£*  / 


>*  * 'rK^QC^ 9W o./u ) 

. -tU] a~[ v (U)-M,(k,)  _ J 

• <**&/£) -*A)  ~f^t]^*[-^H)4  utj 


^ ~ C f>  m O 


(39) 
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hote  that  in  equations  (36)  and  (39)  the  first  double 
intensity  of  the  incident  vase  and  the  second  double  suMatlon  is  the 
intensity  of  the  scattered  vase.  The  last  two  double  aunnations  represent 
interaction  between  the  incident  and  scattered  fields.  The  pretence  of 
Interaction  terns  in  the  expressions  for  the  intensity  components  is  due 
to  the  fact  that  the  total  pressure  and  component  particle  velocities  each 
involve  two  terns.  Upon  aultlplication,  four  products  are  obtained,  two 
of  which  are  cross  product  terns  that  represent  interaction  between  the 
incident  and  scattered  fields.  The  contribution  of  these  terns  to  the 
total  intensity  field  is  significant  when  intensity  neasurenents  are  nade 
close  to  the  bubble.  It  should  be  noted  that  the  total  pressure  and  com- 
ponent particle  velocities  do  not  involve  cross  product  terns  since  these 
can  be  expressed  as  suns  of  the  incident  and  scattered  fields.  In  order 
to  find  the  average  total  radial  and  angular  intensities,  the  tine  averages 
of  equations  (38)  and  (39)  nust  be  date  rained. 

The  average  total  radial  intensity  is  found  to  be 
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The  average  total  angular  intensity  ie  eiven  by 

^ ~ 2pc  ^ 2r0j2z  % (o»  o)fta«Q) ^ 


c*+-[-  JM+  £>' 


w] 


/r  / 


V /-Z  +,)% ft*e)P/(c^ $)A (Ir)^ (% )Mz^ (U) 

•z^w  otp(ka) OX  [*m  (&*)-<=*»,  (kr\  - -^r  r»  — cy  ( l&)  -A^  (ir)-f- 


+£z  1 2 2( ^ *7>  *)  i?mr  ^ o)s  q,  )jt  ar ) 

f m*  C ^ ~ o ' a 

&*)<++[-%*  -Off  (U) 
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i 


+ zZjtZ"  Z )a^  *ju  :> 

I m~0  - O • 

•/^(kr)c *+-i‘X~>ft*)-'*rn  (t:r)  - j»>4-  j 


(1*1) 


The  quantity  is  the  average  intensity  of  the  incident  plane 

wave.  Ifc.g.s.  units  are  used,  the  intensity  is  eeasured  in  ergs  per 
second  per  square  centimeter. 

The  case  of  very  long  wavelengths  is  of  particular  interest. 

The  expressions  for  the  intensities  can  be  simplified  slightly  by  making 
use  of  well  known  asymptotic  formulae.  The  formulae  needed  are  the 
following: 


¥ 

■j*,  [ko) 
n0  (ka. ) 
h„(kA 


(U) 


1-3-5  (2™  + >) 

1 

TaT 

Zr*' — 

Substituting  these  results  into  equation  (11) 
•/ttvu  <*c  (ko~  ) — * (k&0 

, 2r»+  I 


tv  — C 
f»>  0 
nn  ■=?  0 
M > 0 


ka.-*o 

kfr-+o 

ka.  —a 

I 


(la) 


fn  = O 


•U^.«„(ka)—+  - ; — nT 

l /■  3-5  — fort- OIL  !‘3-5  • * • • (2*>- OJ 


( ) 


M >0 


(1*3) 


for  small  value e of  ka-  . Since  the  arguments  of  the  tangents  are  small 


(ko)  (ka. ) - (ko~) 


*Y>  — O 


**m  (Ilk)  ^ A*^cx„(k*)  **  ~ ^ ' 

[l-S-5  /2*n+d]L  /■3-5--(2*i-0]  5 

rr~!  > O 

By  making  use  of  the  above  approximations,  the  terms  in  equations 
(35)  and  (36)  involving  cx^  f J;aS)  can  be  computed  with  ease  for  the  case 
of  long  wavelengths. 

The  average  total  intensity  of  the  entire  sound  field  is  given  by 


fT-lX  + , 


where  th  and  tp  are  the  unit  vectors  in  the  direction  of  increasing  h 
and  & . The  magnitude  and  phase  of  1CT  are 


\%\  = fx;  +rt 


0*7) 


In  the  above  equations  ITrl  is  the  magnitude  of  the  average  total  inten- 
sity at  the  point  £ is  the  angle  between  the  average  total  lntan- 


sivy  and  the  unit  vector  t * • The  angular  deviation  between  the  direction 


of  the  average  total  intensity  produced  by  a plans  wave  incident  upon  a 
spherical  bubble  and  the  average  intensity  of  the  same  plane  wav#  with 
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no  bubble  present  Is  given  by 

A * & + $ * (1,8) 

The  quantity  A will  be  celled  the  Intensity  Deviation  Angle. 
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III.  RESULTS 

The  ratio  of  average  scattered  intensity  to  average  incident 
intensity  has  been  plotted  as  a function  of  angle  for  different  val.uea 
of  /nfi/A.  and  t~/o- . Two  seta  of  curves  have  been  drawn  and  are  presented 
in  Figures  2 and  3.  In  the  first  set  of  curves  2rr */,K  is  held  constant 
and  f/d-  is  varied;  in  the  second  set  of  curves  f"/<L  is  held  constant  and 
2rro/K  is  varied.  The  following  sets  of  values  have  been  taken  for  pur- 
poses of  calculation: 

Figure  2:  2rfa/A-  0.10,  . . } x.  » 3/2,  2,  1* 

Figure  3t  r/a  • 2.00,  2rra/A.»  0.0$,  0.01,  0.02 
For  the  case  in  which  h/xla  fixed,  there  seens  to  be  little  change  in 
the  ratio  of  scattered  intensities  fbr  the  values  of  2rxfK  considered. 
Note  that  for  a wavelength  of  100  feet  and  a bubble  dianeter  of  3 feet, 
the  value  of  is  approximately  0.10.  For  this  value  the  ratio  of 

intensities  1.5  feet  away  froa  the  surface  of  the  bubble  in  the  backward 
direction  along  the  polar  axis  is  0.309. 

The  average  total  sound  intensity  of  the  entire  sound  field  is 
equal  to  the  vector  sun  of  the  average  total  radial  and  average  total 
angular  intensity  covenants.  Unfortunately,  the  nunerical  calculation 
of  these  intensity  components  is  lengthy  and  tedious  since  each  component 
is  equal  to  the  sun  of  four  double  suasaations  in  which  the  indices  run 
froa  zero  to  infinity.  For  co^mtational  purposes,  however,  only  a 
finite  nuafcer  of  these  terns  need  to  be  taken  because  terns  involving 
large  values  of  the  indices  are  ertrenely  naall.  It  is  found  that  for 


I 
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small  nlws  of  the  parameters  2rra//.  and  ''/^sufficient  accuracy  la 
obtained  by  taking  the  largest  ralue  of  each  index  to  Be  two.  Each 
Intensity  component  can,  therefore,  be  expressed  as  the  sua  of  thirty-six 
terms  pro  Tided  calculations  era  Halted  to  the  case  of  long  wavelengths 
and  distances  close  to  the  bubble.  As  the  Talus  of  either  2n*/K  or  /-/*«* 
is  increased,  It  is  found  that  acre  and  acre  terns  are  needed  to  obtain 
the  correct  values  of  the  Intensity  components.  Numerical  work  become 
prohlbltiTs  when  the  wars  length  is  of  the  earn  ordsr  of  magnitude  as  the 
bubble  radiua. 

In  order  to  obtain  torn  Idas  aa  to  the  influence  of  a spherical 
bubble  on  the  a re rage  total  intensity,  calculations  hare  been  performed 
for  ralue s of  2rr&/\*  0.1  and  r/*-«  2,  U,  10.  Them  ralue s are  repre- 
sents tire  of  a small  bubble  in  a sound  field  of  long  ware  length  and  corre- 
spond approxlmetely  to  e bubble  diameter  of  3 feet  end  a ware  length  of 
100  feet.  Figure  U is  e plot  of  the  Intensity  Derietlon  Angle  as  a 
function  of  the  polar  angle  O for  2^/A.  » 0.1  and  r/*-"  2.  Figure  5 
is  e polar  plot  of  the  ratio  of  the  are  rags  total  sound  intensity  to  the 
a re  rage  intensity  of  the  incident  plane  ware  for  2tto./K  ■ 0.1  and  • 2. 
Note  that  the  arrows  shown  on  the  polar  diagram  indicate  the  direction 
of  the  total  sound  intensity  rector  with  respect  to  the  polar  axis  at 
the  point  (2<y,  9 ).  The  ralue  or  irTl/rc  1.5  feet  away  from  the  surface 
of  the  bubble  in  the  backward  direction  along  the  polar  axis  is  0.185. 

Figure  6 is  e polar  plot  of  the  ratio  of  the  are rage  total  sound  intensity 
to  the  a re  rage  intensity  of  the  Incident  plane  ware  for  2**-/Km  0.1  end 
» 2,  U,  10,  oo  . It  is  erldent  that  for  large  rains s of  t/fi-  the  pres- 
ence of  the  bubble  has  little  influence  on  the  are  rage  total  sound  intensity. 
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